LEPTON - PHOTON 2007 



Dark Matter Decaying Now 

Jose A. R. Cembranos, Jonathan L. Feng, and Louis E. Strigari 

Department of Physics and Astronomy, University of California, Irvine, CA 92697, USA 



The instability of dark matter may produce visible signals in the spectrum of cosmic gamma-rays. We consider 
this possibility in frameworks with additional spatial dimensions and supersymmetry. Examples of particles include 
superweakly-interacting massive particles such as gravitinos in supersymmetry models, the lightest Kaluza-Klein (KK) 
^ ^ state in models with universal extra dimensions, and weakly-interacting massive particles such as branons in flexible 

■ brane-worlds. 

^5 Introduction 

< 

■ Though we know that approximately one fourth of the energy density of the Universe is made of non standard 
particles, the microscopic nature of such particles remains as one of the major mysteries in science. This fundamental 
piece of standard cosmology, known as dark matter (DM), is usually assumed to be in the form of stable, massive, 
collisionless and non-self-interacting particles. Weakly-interacting massive particles (WIMPs) that freeze-out with the 



2 ' right thermal abundance in the early Universe are well-motivated dark matter candidates with masses and interaction 
c/2 cross sections of the order of the weak scale. WIMPs emerge naturally from different particle physics scenarios such 



as the lightest supersymmetric particle (LSP) in R-parity conserving supersymmetry (SUSY) models [H, the lightest 
Kaluza-Klein (KK) excitation (LKP) in models with universal extra dimensions (UED) 2], or branons in flexible 
brane-worlds [3|, 



, Other well established candidates are superweakly-interacting massive particles (superWIMPs), which are born 
from unstable WIMPs with typical lifetimes r ~ 1 s - 10* s, and naturally inherit the abundance of the WIMPs from 



■ which they are produced 0, . Examples of superWIMPs include non-thermally produced weak-scale gravitinos 0, 
SQ, SS], axinos [lo| . and quintessinos 11 1 in supersymmetry; and Kaluza-Klein graviton and axion states in models 



od 

■ with universal extra dimensions [12 1 . 

The experimental search for DM depends on its nature, and in general, a strong statement would need the interplay 
• • of collider experiments Il5| and astrophysical observations. These latter observations are typically classified as 
, direct or indirect searches (see [13[ however, for different alternatives). Additionally, elastic scattering of DM particles 
^ ' from nuclei should lead directly to observable nuclear recoil signatures. 

' In this note we point out that present DM particles can be unstable and have associated lifetimes of the order of 
the age of the Universe or longer. In this case, new DM signatures can be explored, such as anomalies in cosmic 
gamma-rays. Such DM particles can be found in the above mentioned scenarios with supersymmety and extra 
dimensions, not only in form of superWIMPs, but also in form of WIMPs, for example in the case of branons. 



Gamma-ray background 

For highly degenerate particles, the decays can take place very late to soft photons. The photon spectrum is not 
thermalized and produces bumps in the diffuse photon spectrum that may be observable. The present differential 
flux of photons from a general decay is given by 

_cn^ f*° N{t)dT^ de^ 
dEry 47r Jo Vq de-y dEj 
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Figure 1: Contributions from extragalactic dark matter decays along with data for the the cosmic gamma background observed 
by COMPTEL. The curves are for decay in mUED with hfetime r = lO^to and 771^1 = SOOGeV (solid) and, B decay 
in SUSY with lifetimes t — 5 x lO^to and frig — 80GeV (dashed). We have assumed that these particles make up all of 
non-baryonic dark matter {Q,nbdm ~ 0.2) and smeared all spectra with energy resolution AE/E — 10%, characteristic of 
COMPTEL [H. 



where n-y is the number of photons produced in a single decay, to — 4.3 x 10^^ s. is the age of the universe, 
N(t) = iV'"e~*/'^, where TV'" is the initial number of decaying particles, Vq is the present volume of the universe, and 
the relation between the produced energy e-y, and observed energy is given by the scale factor of the Universe: 

de-y 



(1 + z) = a ^. In the case of a two body decay, the produced photons have a determined energy given by the 



mass splitting between the decaying particle and the produced one: = AM. It means that 



[i,Ei and 



dE-, 



_ c Arine-P(iJT/eT)/^ / 1 \ 



(2) 



where 



P(a) 



2 (In [\/niF + v^M + f^Afl^j - In [\/n^]) 



(3) 



is given by the equation: 



da 
'dt 



= aH~Hr 



M 



+ n^a^ = Q{a) , 



(4) 



The flux has a maximum at 



M 



(5) 



where U{x) = £+i^v|iH. Here we are neglecting the radiation content, ^ 0, and the curvature: fc ^ [16 
For this calculations, we use a flat cosmological model with il„i = 0.3, J^a — 0.7, and h — 0.7. An interesting side 
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point of this analysis is the sensitivity of the signal to cosmology. The gamma ray background may provide novel 
constraints on the above parameters and dark energy properties. 

We can gain further insight into the spectrum shape in Fig. [T] by considering asymptotic limits. If the lifetime 
is much shorter than the age of the Universe, HqT^Qa « 1, the flux grows as d^/dEj cx E^^^ until it reaches its 
maximal value at: 



jTimax 
^7 



(6) 



From this energy, the flux is suppressed exponentially due to the decreasing number of decaying particles 5]. On 
the other hand, if the lifetime is much longer than the age of the Universe, H^T^Vlf^ >> 1, the flux only grows as 
d^/dE^ oc E^/"^ when the dark matter decays in the matter dominated epoch. For decays in the vacuum dominated 
Universe, the flux decreases as d^/dE^ cx E~^, reaching its maximum at the transition between these two regimes: 



^7 



(7) 



Thus at first order the maximum does not depend on r [l6| . 

In figured! we show the generic region of parameter space Am-r that is probed by the diffuse photon background. 
We use the limit on the diffuse photon background from ^ keV-100 GeV as determined in The solid lines label 
the points where the peak of the photon spectrum from decays exactly matches the observed photon background. Of 
course, for energy regimes where the photon background is well-resolved, these lines provided a conservative upper 
limit to the amount of photon flux allowed from dark matter decays. For a fixed decaying DM (DDM) mass scale, 
all points above these lines are roughly inconsistent with the photon background, as the peak of the decay photon 
spectrum lies above the observational limits. For lifetimes t ^ to, the peak of the energy spectrum is very similar 
to the mass difference. In this limit the relationship between the observed peak and the mass difference is given by 
Eq. ([7]) . The turnover for lifetimes ^ to reflects the fact that the photons of an emitted energy redshift to appear at 
different observed energies. For example, the decay of an 80 GeV particle with lifetime r ~ 10^^ sec. {z ~ 4000) and 
Am ~ MeV produces a spectrum with peak ^ keV. 

The shaded band in figure [2] shows the region where the energy of the detected photons fall in the unaccounted for 
range of 1-10 MeV. The diagonal line shows the relation for the model ^ + 7, though similar relations hold 
for the other gravitational suppressed two-body decays. Importantly, these relations are independent of the mass 
scale for these high degeneracies. The region of parameter space defined by t ~ 10^"'^ sec and Am ~ MeV defines the 
preferred mUED and SUSY models that can account for the cosmic gamma-ray background in the 1-10 MeV regime. 



Gamma-ray lines from the center of the Galaxy 

If the DDM comprises the majority of the dark matter today, they will decay inside our local halo producing 
photons, positrons, electrons and neutrinos. The fluxes of all these particles will reach the Earth providing the 
possibility of a directional dependent DM signal. If this DDM produce photons directly from two-body decays, 
these photons will appear as a line flux at energy exactly determined by the mass splitting, with no redshifting. To 
determine the Galactic signal, we define as the angle between the direction of the Galactic center and the line of 
observation, and / as the distance from the Sun to any point in the halo. In a linc-of-sight direction determined by 
the differential 7-ray fiux from decays of dark matter is 



d$ dN^ 1 



!2 . e-*"/ 



dfldE-y dE^ 47rTma; 



pvmm- (8) 



The radial distance r is measured from the Galactic center, and is related to / by = + Dq — 2I?qIcos\E', where 
Dq = 8.5 kpc is the distance from the Sun to the center of the Galaxy. The distance from the Sun to any the edge 



XXIII International Symposium on Lepton and Photon Interactions at High Energy, 
Aug 13-18, 2007, Daegu, Korea 




log [lifetime (sec)] 

Figure 2: The general mass difference-lifetime parameter space for decaying dark matter scenarios. The solid lines reflect 
the points where the peak of the spectrum from decays exactly matches the observed photon spectrum. The shaded region 
corresponds to a detected energy in the unaccounted for range 1-10 MeV. The diagonal line shows the relation for the LKP 
model in mUED 



of the halo in the direction 9 is Imax = Dq cos 6 + y ~ -C*© sinfl. Each decay produces one photon with energy 
given by the mass difference, Am, so the photon spectrum is dN^/dE^ — 5{E^ — Am). 

The photon flux is maximized in the direction of the Galactic center, and must be averaged over the solid angle 
of the detector. For detectors with sensitivities to energies in the tens of MeV regime, the solid angles are typically 
of order Afi = 27r(l — cos^*) ~ 10~^. One can model the dark halo of the Milky Way with the results of !l8|. We 
have found that the effects of dark matter clumps in the halo are negligible. Indeed, for realistic distributions of the 
clumps, the flux is only increased by ~ 1% for all lifetimes. 

Positrons may also be produced in the decays. If positrons are produced with energies smaller than 3MeV, 
they will be effectively stopped before annihilating with electrons, losing their kinetic energy through coUisional 
ionization or excitation in neutral Hydrogen and by interaction with plasma waves in ionized interstellar medium. 
In such a case, they produce a narrow photon line at ~ nif, ~ 511keV if they annihilate directly into 2 photons 
or through parapositronium formation. Annihilation in such a state happens 25% of the time that positronium is 
formed. The other 75% of the time, the annihilation takes place in the orthopositronium state, yielding a photon 
continuum with < 511 keV. Consequently, the total number of 511 keV photons produced per unit time is given 
by dn^ii/dt = 2(1 — 3p/4)n£i£ij\/rg+ = 2(1 — ip/ A)pbbm^ e+ 1 ^ where we are supposing that annihilation takes 
place through positronium formation a fraction p of the time. 

The predicted distribution of the 511 line for any particular model of the Milky Way dark halo can be computed 
as the integral along the line of sight, as a function of galactic longitude I and latitude 6, of the emissivity dn^u/dt: 

dE-ydn ATTjf) dt ^ ' 

where in this case, we are parameterizing the halo radius as: r{s,b,l) — Dq + + 2Dqs cos b cos I. The morphology 
of the emission depends on the Milky Way halo profile (in addition to the precise distribution of baryons and their 
chemical composition, which determine the positron propagation). Observations are just sensitive to the inner region 
of the halo, and demand a cuspy halo. The intensity of the emission is so low in outermost regions that it is difhcult 
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to discriminate from the instrumental background. Interestingly, the splittings and lifetimes that can account for the 
cosmic gamma-ray background at MeV energies are also appropriate to explain the 511 line. An example of a DM 
candidate that can provide such phenomenology is the branon, due to its universal coupling with all standard model 
particles [19 1. 



Conclusions 



The main idea of this note is that decaying dark matter can provide indirect detection signals at scales much smaller 
than their masses. We have analyzed this idea for the extragalactic photon background observed by COMPTEL [2^ 
and the photon lines coming from the Galactic center observed by INTEGRAL 2l|. Below ^ 1 MeV and above 
^ 10 MeV, the extragalactic photon fluxes are well-modeled by the diffuse cosmological emission of unresolved Active 
Galactic Nuclei |23] , but in the entire regime of ^ 1 — 5 MeV, no population of known astrophysical sources are able to 



account for the observed cosmic gamma-ray background 2^ 24 1 • An additional intriguing source of non-relativistic 
positrons from the center of our Galaxy has been determined recently by the INTEGRAL satellite. 

We have shown that DDM provides a viable explanation that deserves further investigation. As we have seen, 
these late decays at MeV scales are natural for gravitationally suppressed superWIMPs, and also for TeV WIMPs 
such as branons in flexible brane-worlds. 
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